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Abstract 

We have developed a numerical model for the tempo- 
ral evolution of particle and photon spectra resulting 
from nonthermal processes at the shock fronts formed 
in merging clusters of galaxies. Fermi acceleration 
is approximated by injecting power-law distributions 
of particles during a merger event, subject to con- 
straints on maximum particle energies. We consider 
synchrotron, bremsstrahlung, Compton, and Coulomb 
processes for the electrons, nuclear, photomeson, and 
Coulomb processes for the protons, and knock-on elec- 
tron production during the merging process. The 
broadband radio through 7-ray emission radiated by 
nonthermal protons and primary and secondary elec- 
trons is calculated both during and after the merger 
event. To test the ability of the computer model to 
accurately calculate the nonthermal emission expected 
from a cluster merger event, we apply the model to the 
Coma cluster of galaxies, and show that the centrally 
located radio emission and the Hard X-ray excess ob- 
served at 40-80 keV is well fit by our model. If our 
model is correct, then the Coma cluster will be sig- 
nificantly detected with GLAST and ground-based air 
Cherenkov telescopes. 

1 Introduction 

Merger events between clusters of galaxies are ex- 
tremely energetic events. The formation history of a 
cluster of galaxies will include several merger events 



( Gabici & Blasil 2003 ). With typical masses for a 
rich cluster ~1O 15 M , dynamical estimates of the en- 
ergy deposited in to the i nternal s tructur e of a clus- 
ter is ~ 10 63 -10 64 ergs dSarazinl |2004) . Approxi- 
mately 5% of the available energy is assumed to ac- 
celerate particles from the thermal pool to create a 
distribution of nonthermal particles. The first-order 
Fermi process is capable of accelerating particles to 
~ 10 19 eV by the shocks that will form at the inter- 
action boundary of two m erging cluster of galaxies 
I Berrington & Dermeii 12003). 

Optical and X-ray studies estimate that ~ 30-40% of 
rich cl usters show signs of a current or recent merger 
event dForman et al. . 1981 : Beers. Geller. & Hu chra. 



1982). Structure formation calculations estimate that 



a cluster will s ee several mergers t hough out its for- 



mation history dLacev & ColeL 119931) . indicating that 
most clusters contain a population of highly energetic 
nonthermal particles. 

Diffuse cluster radio emission with no associated com- 
pact counterpart is observed in an ever increasing 
number of galaxy clusters, providing evidence for the 
prevalence of diffuse, highl y energetic nonthermal par- 
ticles in galaxy clusters. Jg iovan nini. To rdi. & FerettiL 
19991 Kempner & Sarazin, 2001) The diffuse radio 
emission is classified in to two categories. Diffuse 
emission found in the cluster center that mimics the 
thermal bremsstrahlung emission with random polar- 
ization is known as a radio halo. Diffuse cluster emis- 
sion located on the cluster periphery characterized by 
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irregular shapes and strongly polarized light is known 
as a radio relic. These features are preferentially seen 
in clusters with current or recent cluster merger events, 
and are thought to be the observational consequences 
of nonthermal particles accelerated by shocks formed 
in a cluster merger event. 

2 Models 



We present results of a 



( Berrington & Dermer 
late the time-dependent evolution 



numerical model 
2003) designed to calcu- 
of nonthermal 



particle distribution functions evolving through ra- 
diative losses. The electrons and are accelerated by 
the first-order Fermi process and at the shock fronts 
formed in merging clusters of galaxies, and the result- 
ing photon radiation is calculated. Particle injection 
functions Q e ,p f° r electrons ("e") and protons ("p") 
are assumed to be described by power-law momentum 
spectra. In terms of kinetic energy K ep , the injection 
functions are given by 



Q° [K etP (K eiP + 2m eiP c 2 )l 



s+l . 

~~ and 



(K e , p + m e)P c 2 ) exp 



K, 



e,p 



(1) 



where s is the injection index and K m3iX is the maxi- 
mum particle energy determined by three conditions: 
the available time to accelerate to a given energy since 
the formation of the cluster merger event; the require- 
ment that the particle Larmor radius is smaller than 
the size scale of the system; and the condition that the 
energy-gain rate through first-order Fermi acceleration 
is larger than the energy-loss rate due to synchrotron 
and Compton processes. The constant Q ep normalizes 
the injected particle function, and is determined by 

ax 



dK ep K ep Q ep (K ep ,t) 



Ve,p 



A r]n e m p vl{nicM) , 



(2) 



where (nicivi} is the number density of the intra-cluster 
medium (ICM) averaged over the area A of the shock 
front, r/y c is an enhancement factor to account for the 
ions heavier hydrogen, 7] eiP is an efficiency factor taken 
to be 5% unless otherwise noted, and m p is the mass 
of a proton. The minimum kinetic energy K mm is held 
constant at 10 keV. 

The model calculates the forward [v\] and reverse [^2] 
shock speed from the gravitational infall velocity v g 



of a merging cluster. The trajectory of the smaller 
merging cluster is approximated by a point mass of to- 
tal mass M2 that falls onto a dominant cluster of to- 
tal mass Mi whose density profile is described by an 
isothermal beta model. The velocity v of the shock 
fluid is calculated by solving the equation 



/xi ni 

V-2 n 2 



1 + 3M± 

1 +2,Mn 



(3) 



where ni and n 2 are the number densities in the dom- 
inant and merging cluster, respectively. The mean 
atomic mass in the dominant cluster and the merging 
cluster are given by n\ and H2, respectively. Both mean 
atomic masses are set equal to 0.6m p . 

The Mach speeds of the forward [.Mi] and reverse 
[M.2] shocks are calculated by 



Mi 




2v n 



(>2 



1 + 1/1 + 



4 k 



(4) 



where c\ is the sound speed in the dominant cluster, 
and C2 is the sound speed in the merging cluster. The 
Mach number is defined to be At 1,2 = ^1,2/01,2 for the 
forward and reverse shock, respectively. Equation [5] 
and equation|4]are derived from the shock jump condi- 
tions, by equating the energy densities of the forward- 
and reverse-shocked fluids at the contact discontinuity. 
Compression ratios C\ (forward) and C2 (reverse) are 
calculated from the equation 



C 



1,2 



3(1 



Mil) 



(5) 



The time-dependent particle spectrum N(K,t) is de- 
termined from solving the Fokker-Planck equation in 
energy space, given by 



dN(K,t) 1 d 2 



Of 



2dK 2 
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(6) 
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Figure 1: Light curves at various observing frequencies 
produced by a shock formed in a merger between 10 14 M© 
and 10 15 M Q clusters that begins at Zi — 0.3 and is evolved 
to the present epoch (t = 3.42 Gyr). All light curves are 
for a magnetic field strength of B = l.OfiG unless other- 
wise noted. Radio light curves in Jansky units are given at 
15 MHz, 30 MHz, 74 MHz, 110 MHz, and 1.4 GHz. The 
15 MHz light curve is also calculated with a magnetic field 
strength of B = 0.1/iG (dotted curve). 

Here K(K, t) is the total kinetic-energy loss rate found 
by the sum of Coulomb, synchrotron, bremsstrahlung 
and Compton processes for electrons, and Coulomb 
processes for protons. In addition, the protons expe- 
rience catastrophic losses from proton-proton (i = pp) 
collisions, proton-7 (i = p^y) collisions, and diffusive 
escape (i = d) on the timescale Ti(K,t). Secondary 
electrons are calculated and added to the primary elec- 
tron distribution function and are subject to the same 
energy losses as the primary electrons. 

3 Results 

Light curves at various observing frequencies are 
shown in Figures [2 & El using the parameters for the 
standard case of a merger event between a dominant 
cluster of mass Mi = 1O 15 M , and a merging cluster 
of mass M2 = 10 14 M Q with a magnetic field strength 
B = l.OfiG, and a beginning redshift of Z{ = 0.3. 
The light curves of the nonthermal radiation exhibit a 
common behavior independent of frequency. At early 
times, the spectral power rises rapidly as the clusters 
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Figure 2: Same as Figure ^except light curves in energy 
flux units are given at 40-80 keV, >100 MeV, >1 GeV and 
>100 GeV. Curves are labeled in the diagram. 



merge. The peak emission occurs when the centers 
of mass of the two clusters pass at t co n, after which 
the emission exhibits a slow decay and approaches a 
plateau at times t > i acc . The time i acc .i(2) is defined 
as the time at which particle injection has terminated at 
the forward shock (1) or reverse shock (2). The rate of 
decay of the emission increases with radio frequency 
due to the stronger cooling of the higher energy elec- 
trons, so that the decay is slowest at lower frequen- 
cies. Synchrotron emission from secondary electrons 
forms the late-time plateaus at radio energies. This be- 
havior is also apparent for the hard X-ray emission, 
although it is formed by primary bremsstrahlung and 
both primary and secondary Compton radiation at late 
times. At 7-ray energies, the 7r°-decay emission forms 
a plateau of emission that dominates soon after i aC c- 

Calculations of the hardest particle injection spectral 
index s m - m formed in cluster merger shocks are shown 
in Figure |3] as a function of the larger mass Mi of the 
two clusters, with a constant subcluster mass M2 = 
10 14 M Q . We also assume that the onset of the merger 
begins at redshift Zi = 0.1; softer injection indices 
are obtained for mergers at larger values of Zj, because 
of the smaller maximum separations for merger events 
occurring at higher redshift. The mean cluster masses 
are smaller at higher redshift. These lower mass clus- 
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Figure 3: Calculations of the hardest particle injection 
spectral indices s m i n formed in cluster merger shocks as 
a function of the larger mass M\ of the two clusters, for 
various values of r c and (3. The values of r c are given in 
Mpc. The minimum spectral index for the forward and re- 
verse shock is shown in the figure. 



Figure 4: Comparison of radio observations of the Coma 
cluster with cluster merger shock model. The solid curve 
is the total expected radio emission, the dashed curve is 
the contribution to the total emission from primary elec- 
trons, and the dotted curve from secondary electrons and 
positrons. The solid circles are observational d ata points 
taken from lThierbach. Klein. & W ielebinskil j2003l) . 



ters have a lower virial temperature. If this is consid- 
ered in the calculation of s m i n , then it is possible that 
stronger shocks may be seen at higher redshifts. We 
calculate s min for various values of r c and f3. The val- 
ues (r c ,/3) = (0.05, 0.8), (0.05, 0.45), (0.5, 0.8), and 
(0.5, 0.45) roughly correspond to the extrema in the 
range of these parameters m easured for 45 X-r ay clus- 
ters observed with ROSAT dWu & XueL l2000h . Also 
shown are values of s m - m for (r c , j3) = (0.179, 0.619), 
which are the average values of these parameters for 
the 45 X-ray clusters, and (r c ,(3) = (0.25,0.75), 
which are the standard parameters used in the calcu- 
lations. It should be noted that as the dominant cluster 
approaches the mass of the merging cluster (1O 14 M ), 
the spectral indices for the forward and reverse shock 
will be identical if the matter profiles are identical. 
However, in Figure[3]the values of (r c , j3) are held con- 
stant, and therefore; the matter profiles differ. 

3.1 Modeling the Coma Cluster 



Both opt ical dColless & Dunrl ll99ot biviano et al 
1 1994 Edwards et all l2002h and X-ray 
Jvikhlinin. Forman. & Jonesl 1 19971: lArnaud et al 



2001) observations indicate that the dynamics of the 



Coma c luster is well-described b y a three-body merger 
model. Colless & Dunnl ( 199ot) were the first to find 
substructure in the central region of Coma that is 
consistent with a recent merger event near the collision 
time t co \\, defined to be the time when the centers 
of mass of the two clusters pass through each other. 
Total mass estimates of the dominant cluster M\ is 
estimated to be 0.8 x 10 15 M^. X -ray observations 
Jvikhlinin. Forman. & Jonesl 1997 ) estimate the total 
mass of the merging cluster M2 from the gas striped 
in the merging process to be ~0.1 x 10 15 M Q . This 
assumes a gas fraction of 5-10%. 

The ICM is well described by an isothermal beta 
model with core radius r c = 0.257 Mpc, central elec- 
tron density p e Q = 3.82 x 10~ 3 cm~ 3 , central pro- 
ton density po = 7 A3 x 10~ 27 gm cm -3 , power- 
law slope = 0.705, and a mean gas tempera- 
ture (T x ) = 8.21 keV (Mohr. Mathie sen. & Evrard . 
1999b . The assumed magnetic field strength is B = 
0.22/xG, and an efficiency factor r/ ep = 1%. Be- 
cause we can never know the true gas distribution 
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Figure 5: Comparison of X-ray observations of the Coma 
cluster with the cluster merger shock model. Points are la- 
beled in the diagram. The thermal bremsstrahlung emis- 
sion (dashed curve) is the calculated thermal bremsstrahlung 
for the parameters described in the text. The OSSE data 
points are 2a upper limits. The solid curve is the sum of 
the non-thermal (dotted curve) and the thermal 8.21 keV 
bremsstrahlung emission. 



in the merging cluster, we approximate the ICM of 
the merging cluster by an isothermal beta model with 
core radius r c = 0.150 Mpc, central electron den- 
sity p e0 = 1.0 x 10~ 3 cm" 
sity pq = 1.67 x 10 -27 gm cm" 
= 0.7, and a mean gas temperature calc ulated from 
equation 9 in lBerrington & Dermen 12003). 

The model is evolved to a time t = 0.97 x 10 9 yrs, 
which is just prior to the collision time t co \\ = 1.0 x 
10 9 yrs. The redshift of the cluster at the time corre- 
sponding to the creation of the shock front corresponds 
to Zi = 0.10 for a (f2o,£X\) = (0.3,0.7) cosmol- 
ogy. The observed redshift of the evolved cluster is 
z = 0.0232. 

In Figure |4j we show a comparison of our model 
with the observed radio emission from the radio halo 



central proton den- 
" 3 , power-law slope 



Thierbach, Klein, & Wielebinski (2003). Our models 
favor a primary electron source for the radio emis- 
sion from the radio halo. Despite our models using 
a uniform density profile for the calculation of the sec- 
ondary electron production, the emission from the sec- 



ondary electron is an upper limit to the true secondary 
electron emission. However, the density in the central 
region is roughly constant, so that the uniform density 
assumption is a good approximation to the true sec- 
ondary electron emission. 

In Figure |3J we show a comparison of the calcu- 
lated thermal and non-thermal emission with the ob- 
served Hard X-ray (HXR) emission from the cen- 
tral region of the Coma cluster of galaxies observed 



by the Phos wich Detection System (PDS) on Bep- 

|2004) . The reported 



poSAX ( Fusco-Femian o et al 



non-thermal photon emission is the total integrated 
emission expected within 1.5 Mpc of the cluster cen- 
ter. This corresponds to a field of view of ~ 1?7. 
The PDS has a FWHM field of view of ~ 1?3 which 
corresponds to a linear scale of ~ 2.2 Mpc. In 
addition we also show the OSSE 2a upper limits 
( Rephaeli. Ulmer. & Grubel Il994b . The HXR emis- 
sion observed between 20-80 keV is dominated by 
Compton up scattering of photons off of primary elec- 
trons. 

The BeppoSAX/PDS observations include the merging 
cluster associated with NGC 4839, and will be contam- 
inated with any non-thermal emission associated with 
shocks formed in its merging process. The expected 
non-thermal emission resulting from shocks for a clus- 
ter infalling at a mini mum distance of l-6h§ Mpc 
( Neumann et all 1200 lb will be negligible in compari- 
son with the e mission from the merger obse rved in the 
core of Coma JBerrington & Dermei , 2003 ). 

In Figure |6j we present the predicted 7-ray emission 
from the Coma cluster of galaxies. The ob servational 
limits are taken from Iweekes et al ( 2002b . As seen 



from the figure, the predicted 7-ray emission falls 
comfortably below the predicted upper limits for the 
EGRET observations. Our model predicts that the 
space-based observatory GLAST will detect the non- 
thermal 7-rays at high significance. Furthermore, we 
predict that both VERITAS and HESS will have strong 
rj 5(7 detections in 50 hours of observations. 

4 Summary and Conclusions 

We describe a computer model designed to calculate 
the nonthermal particle distributions and photon spec- 
tra resulting from nonthermal processes produced by 
shocks that form between merging clusters of galax- 
ies. Over the lifetime of a cluster merger shock ~10 61 - 
10 62 ergs will be deposited in the energy of a nonther- 
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Figure 6: Predicted 7-ray emission from the Coma clus- 
ter of galaxies from the cluster merger shock model. The 
solid curve is the integrated number of particles per cm 2 
per s 1 greater than E pca k. Observational limits for EGRET, 
GLAST, and VERITAS and HESS are included. The 
EGRET and predicted GLAST observational limits are 1 
year integrations. The quoted VERITAS an d HESS limits 
are 50 hour, 5<r limits ( Weekes et al., 120021) . The vertical 
axis represents the number of particles greater than E pea k in 
units of cm~ 2 s _1 . 



mal particle population. From the calculated peak lu- 
minosities shown in Figures [2 & El we estimate that 
the nonthermal emission from merger shocks will be 
detected at radio frequencies with a Long Wavelength 
Array out to a distance of ~2000 Mpc at 15 MHz and 
~ 700 Mpc at 120 MHz. At 7-ray energies, we es- 
timate the distance threshold for detection of cluster 
merger shocks to be ~ 200 Mpc for the space-based 
observatory GLAST. However, each cluster must be 
considered on a case by case basis to determine its ob- 
servability. 

It was suggested that Galaxy clusters are a dominant 
contributor to th e diffuse extraga l actic 7-ray back- 
ground (DEGB) (Loeb & WaxmanL l2000h . Our mod- 
els do not support this claim. Shocks formed in the 
merger process are weak resulting in a spectral index 
that is softer than the observed 2.30 ± 0.03 index for 
the DEGB. The superposition of softer spectra on the 
harder spectra observed for the most massive merger 
events will result in a concave spectrum at lower ener- 
gies that is not observed. We estimate that the contri- 
bution to the DEGB by cluster merger shocks will not 
exceed ~1-10%. 

We applied this model to the Coma cluster of galax- 
ies. We show that the radio emission from the ra- 
dio halo Coma C is well fit by the cluster merger 
model. The calculated nonthermal X-ray emission 
also fits the observed HXR emission observed at 40- 
80 keV by BeppoSAX/PDS. The 7-ray emission ex- 
pected from this model is also calculated and is shown 
to fall below the observational limits for EGRET as re- 



ported by Re imer et alJ d2003h . but should be strongly 
detected by the space-based observatory GLAST and 
the ground-based air Cherenkov telescopes such as 
VERITAS and HESS. Even though other acceleration 
mechanisms or point sources could produce nonther- 
mal emission in the core of Coma, our model of clus- 
ter merger shocks account for the entire observed radio 
and HXR emission. 



5 Acknowledgments 

The work of CD is supported by the Office of Naval 
Research and GLAST Science Investigation Grant No. 
DPR-S-1563-Y. 



References 



Arnaud, M., et al. 2001, A&A, 365, L67 

Berrington, R. C. & Dermer, C. D. 2003, ApJ, 594, 709 



X-Ray and Radio Connections http://www.aoc.nrao.edu/events/xraydio Santa Fe NM, 3-6 February 2004 1 



Biviano, A., Durret, R, Gerbal, D., Le Fevre, O., Lobo, 
C, Mazure, A., & Slezak, E. 1996, A&A, 311, 95 

Beers, T. C, Geller, M. J., & Huchra, J. P. 1982, ApJ, 
257, 23 

Colless, M. & Dunn, A. M. 1996, ApJ, 458, 435 
Edwards, S. A., Colless, M., Bridges, T. J., Carter, D., 

Mobasher, B., & Poggianti, B. M. 2002, ApJ, 567, 

178 

Forman, W., Bechtold, J., Blair, W., Giacconi, R., van 
Speybroeck, L., & Jones, C. 1981, ApJL, 243, L133 
Fusco-Femiano, R., Orlandini, M., Brunetti, G., Fer- 
etti, L., Giovannini, G., Grandi, P., & Setti, G. 2004, 
ApJL, 602, L73 
Gabici, S. & Blasi, P. 2003, ApJ, 583, 695 
Giovannini, G., Tordi, M., & Feretti, L. 1999, New As- 
tronomy, 4, 141 
Kempner, J. C. & Sarazin, C. L. 2001, ApJ, 548, 639 
Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 
Loeb, A. & Waxman, E. 2000, Nature, 405, 156 
Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 
517, 627 

Neumann, D. M., et al. 2001, A&A, 365, L74 
Reimer, O., Pohl, M., Sreekumar, P., & Mattox, J. R. 

2003, ApJ, 588, 155 
Rephaeli, Y., Ulmer, M., & Gruber, D. 1994, ApJ, 429, 

554 

Sarazin, C. L. 2004, (these proceedings) 

Thierbach, M., Klein, U., & Wielebinski, R. 2003, 

A&A, 397, 53 
Vikhlinin, A., Forman, W., & Jones, C. 1997, ApJL, 

474, L7 

Weekes, T. C. et al. 2002, Astroparticle Physics, 17, 
221 

Wu, X.-P. & Xue, Y.-J. 2000, ApJ, 542, 578 



